Birefringent break up of Dirac fermions in a square optical lattice 
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We generalize a proposal by S0rensen et al. [Phys. Rev. Lett. 94, 086803 (2005)] for creating an artificial 
magnetic field in a cold atom system on a square optical lattice. This leads us to an effective lattice model 
with tunable spatially periodic modulation of the artificial magnetic field and the hopping amplitude. When 
there is an average flux of half a flux quantum per plaquette the spectrum of low-energy excitations can be 
described by massless Dirac fermions in which the usually doubly degenerate Dirac cones split into cones with 
different "speeds of light" which can be tuned to give a single Dirac cone and a flat band. These gapless 
birefringent Dirac fermions arise because of broken chiral symmetry in the kinetic energy term of the effective 
low energy Hamiltonian. We characterize the effects of various perturbations to the low-energy spectrum, 
including staggered potentials, interactions, and domain wall topological defects. 
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With the discovery of graphene [Q]] and topological insu- 
lators J2l there has been much recent interest in systems in 
which low energy excitations can be described using Dirac 
fermions. A parallel area of interest has been the exploration 
of the possibility of generating artificial magnetic fields for 
cold atoms confined in an optical lattice. Neutral bosonic cold 
atoms cannot couple to a magnetic field directly, so there have 
been numerous proposals 13« 6] of approaches to couple atoms 
to an artificial magnetic field, several of which have been im- 
plemented experimentally 0, 01 ■ 

The problem of the spectrum of quantum particles in a uni- 
form magnetic field on a lattice has the well-known Hofstadter 
spectrum @|. Our modification of the proposal by S0rensen 
et al. lit] leads to an effective Hamiltonian with a tunable 
Hofstadter-like spectrum that arises from the combination of 
hopping and an artificial magnetic field with a non-zero av- 
erage that are both periodically modulated in the x and the y 
directions. The presence of spatial periodicity in the ampli- 
tude as well as the phase of the hopping is the key difference 
between the model we consider here and previous work on 
the spectrum of particles in the presence of magnetic fields 
that are periodic in both the x and y directions fioll . This 
difference facilitates the unusual Dirac-like spectrum that we 
discuss in this Letter. 

In our effective model, when there is an average of half 
a flux quantum per plaquette, and at half-filling, the low en- 
ergy degrees of freedom can be described by a Dirac Hamilto- 
nian with the unusual property that chiral symmetry is broken 
in the kinetic energy rather than via mass terms. This has 
the consequence that the doubly degenerate Dirac cone for 
massless fermions splits into two cones with tunable distinct 
slopes, analagous to a situation in which there are two speeds 
of light for fermionic excitations, similar to birefringence of 
light in crystals such as calcite. We discuss the meaning of 
broken chiral symmetry in our effective model and explore the 
effects of various perturbations, such as staggered potentials, 
domain walls, and interactions between fermions. 

The approach to obtain an artifical magnetic field for cold 
atoms in an optical lattice suggested by S0rensen et al. Jit] 



was presented in the context of the Bose-Hubbard model, but 
ignores interactions and is not specific to bosons. We con- 
sider a model of spinless fermions (corresponding to only one 
available hyperfine state for cold atoms) with Hamiltonian 



H 



(1) 



where c] and c, are fermionic creation and annihilation oper- 
ators respectively at site i, hi = cfc, is the number operator, 
and the notation indicates that we restrict the sum in 
the hopping term to nearest neighbours only. There can be 
no Hubbard-like interaction for spinless fermions, and since 
nearest neighbour interactions in an optical lattice system are 
weak, we postpone our discussion of interactions. 

In Ref. Jfl, two steps are required to generate an artificial 
magnetic field. First, a time-varying quadrupolar potential 
V(t) = V qp sm(uit)xy is applied to the system, and second, 
the hopping is modulated as a function of time. During the 
course of one oscillation of the quadrupolar potential, hop- 
ping in the x direction is turned for a very short period of 
time r <C to = 77 at times t = Tito, where n is an integer, 
and hopping in the y direction is turned on for time r around 
t = (n + h) to. Due to the periodic oscillation in the Hamil- 
tonian, the time evolution operator after m periods may be 
written as U(t = mt Q ) = U(t = t ) m . 




FIG. 1 : Time dependence of the hopping and the quadrupolar poten- 
tial during the course of one period of the quadrupolar potential. 
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Our modification to the proposal in Ref. ]6[] is that when 
hopping is turned on in the a>direction at time t = nto, hop- 
ping is also turned on in the y-direction with an amplitude 
< (3 < 1 relative to the hopping in the x-direction. At time 
t = (n + I J to, hopping is turned on in the y direction, and 
hopping in the .T-direction is turned on with amplitude (3 rela- 
tive to the hopping in the .T-direction as illustrated in Fig.Q] 

The operators for hopping in the x and y directions, 



are T, 



3 ^2 X 



h.c. and T„ 



-x,y -r u.-i'-j respectively, and we may 
write the time evolution operator as 



U(t = mt ) 



e - % ( f * +P f y ) e 2 ™** v e - tt (f>T* + f y ) 

xe -2-rriaxy e -^j-(f :c +pfy) ^ q\ 



where a — V qp /Trhu) and we have set the lattice constant to 
unity. To lowest order in Jt jh we can write this in the form 



[(Jlii; 
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FIG. 2: Spectrum as a function of a when fi = 1 obtained by exact 
diagonalization on a 97x97 site lattice: there is no artificial mag- 
netic field, yet due to the periodic hopping, the spectrum has some 
similarities with the Hofstadter spectrum. 



U = e f 

where the effective Hamiltonian is 

H eS = 



-JoE{[(i+^ m )4 + i^ s + k 

x,y 

+ [(/3 + e 2 "^) £l +hy c x , y + h.c] } , (3) 

with Jo = rJjtQ. A more conventional way to write this 
Hamiltonian is in the form 



(4) 



from which we may identify the amplitude of the hopping: 



tx+l,y 
to 



Joyl + P 2 + 2/3cos(27ray), 
<>x,v+i = Joy/l + P 2 + 2/3 cos(27raa;) , 

and the artificial magnetic field 

2nah { 1 + /3cos(27rax) 



B, = 



1 



1 



(3 2 + 2/3 cos(27raa;) 
1 + j3 cos(27ra?;) 1 
h /8 2 + 2/3 cos(27rm/) J 



(7) 



This field is the sum of a spatially uniform piece with mag- 
nitude and a piece that is spatially periodic in both the 
x and y directions. If j3 = 0, the hopping amplitude is Jo 
and the field is uniform with strength 2v ^ a t corresponding to 
a flux of a0o per plaquette (where <fio is the flux quantum) as 
found in Ref. [6fl and there is a Hofstadter spectrum. If /3 = 1, 
then B z = 0, but the hopping parameters are still spatially 
periodic. At /3 intermediate between and 1, both the hop- 
ping and the magnetic field are spatially periodic in x and y. 
This illustrates the essential difference between the model we 



consider and previous work on quantum particles in a peri- 
odic magnetic field on a lattice - there is spatial periodicity of 
1 1 a in the amplitude of the hopping as well as in the magnetic 
field. For finite /3 the spectrum (illustrated for (3 = 1 in Fig.|2]i 
as a function of a is reminiscent of the Hofstadter spectrum. 

When o = 1/2 there is an average of half a flux quantum 
per plaquette, and the theory is time reversal symmetric as 
fermions cannot detect the sign of the flux 111 ill . The effective 
Hamiltonian Eq. ((3) simplifies to a tight binding model with 
four sites in the unit cell as shown in Fig. [3] a). 





B 


-J_ 


D 


(5) 


J + 


h 


J_ 


(6) 


A 




C 





E 








/ J * 





















FIG. 3: a) Unit cell of tight binding model with hopping parameters 
indicated, b) Dirac cones corresponding to J+ and J_ bands. 

Labelling the four sites in the unit cell as A, B, C, and D, 
and Fourier transforming in space, we may rewrite the effec- 
tive Hamiltonian in the following form: 



H 



= -"W*, (8) 

k 



with 



n k = 2 



/ o 



V 



J + cos k y 
J + cos k x 
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where J± = Jq{1 ± /3) and (in row vector form) ipk 
(cAk,c B k,cck, c Dk ). We find that the dispersion is 



±J± 



and we note that when /3 = 1, J_ = 0, so there will be a 
flat band at e = and a dispersing band associated with J + . 
We can also see that in the vicinity of the points K± ; ± = 
(±f, ±§)> the spectrum is linear 



E q = ±J± 



q 2 y , 



(9) 



where q = k (±-j,±^J, and there are cones with two 
different slopes, corresponding to J± respectively as illus- 
trated in Fig. 0b). When j3 = 0, the two slopes are identical, 
whereas as /? — s- 1, the J_ band becomes flat, and the J+ band 
remains as a cone. Several authors recently considered lattice 
models for cold atoms that are equivalent to the j3 = 1 limit 
of our model, in which there are three bands, one flat, and one 



Dirac like 1121 11311 . When f3 ^ 1 , the underlying Dirac struc- 
ture of the problem is exposed, allowing us to understand this 
unusual dispersion from a symmetry point of view 1 1-411 . 

To start in this direction, we expand around the Dirac points 
and represent the low energy theory (with k measured with 
respect to K) 

H k = 2J [( 7 °7 1 + */?7 3 ) kx + (7°7 2 + */?7 5 ) ky] , (10) 



where we use a non-standard representation of the gamma ma 
trices in which 7 ~ 03 (g) 03, 7 1 = ia 2 ® I2, 7 2 

i/3 J/T. ^ T~ nnA 



i(T 3 (X) cr 2 , 

j" = —io\ Cg) I2, and j b = —jV-f 1 ^ 2 ^ 3 — — i<7 3 g) o\. 
The matrices 7 , 7 1 , 7 2 and j 3 satisfy the Clifford algebra 
7^7" + 7^7^ = 2g^ v with Minkowski metric g^". 

The dimension of the minimal representation of the Clif- 
ford algebra in 2+1 dimensions is 2, allowing for the 2x2 
Pauli matrices as a choice for the 7s. A non-minimal 4x4 
representation as we have used above leads to a freedom in 
the choice of the 7 matrix, i.e. a matrix with (7 ) = 7 4 
that anticommutes with j 1 and r ) 2 . Candidates for 7 are then 
{ 7 °, 7° 7 3 , 7 °7 5 , 7 1 7 2 }- The matrices { 7 °, 7 ° 7 3 , 7 °7 5 } form 
a triplet and r ) lr ) 2 forms a singlet with respect to the SU(2) 
"chiral"-symmetry group with generators {§7 3 , §7 5 , |7 35 } 
(where 7 s5 = 7 3 7 5 ). Each different choice of 70 corresponds 
to a different labelling of the four sites in the unit cell. The el- 
ements of the chiral group generate transformations between 
each labelling. For example, the generator 7 s translates the 
plaquette indices to the labelling of the neighboring lattice 
cell along the y-direction, whilst 7 s translates the plaquette 
indices to the neighbouring cell in the ^-direction. 



(c A \ 

CB 
CC 

\c D J 



CB 
CA 
-CD 

-cc ) 



( c A \ 

CB 
CC 

\c D J 



(c c \ 

CD 
CA 

\c B J 

(11) 



Similarly, 7 35 translates the plaquette one lattice cell along the 
x- and one lattice cell along the y- direction. When (3 = 0, 
the elements of the chiral group are symmetries of Hk- 

When (3 ^ 0, the 7 s and 7 s terms in Hk break the chiral 
symmetry and shifts along either the x- or y-directions do not 
leave Hk invariant. This manifest chiral symmetry breaking 
is inherently different from the conventional notion of spon- 
taneous chiral symmetry breaking in field theoretical models 
which is the signature of mass generation II 1 511 . 

An additional discrete symmetry of the Hk (that arises from 
the hopping structure in Hk) that holds even when /3 ^ is 

r= l( 7 V+7 2 7 5 )-^(7 2 7 3 -7 1 7 5 ), 

which corresponds to a reflection about the diagonal AD in 
the unit cell, with c A ^ c A , cb — > cq, cc —> cb and cd — > 
—cd- The action of T on Hk is to exchange k x and k y . 

Fermion birefringence: As illustrated in Fig. [3] b) the 
dispersion Eq. (O admits massless fermions with two dif- 
ferent "speeds of light" controlled by (3. The eigenvec- 
tors (written as row vectors) for the positive and negative 
energy J + bands are \I>i = ^ (1, — sin9, — cos 9, 0) and 

\&2 = ^75 (lj sin 9, cos 0,0); whilst the eigenvectors for the 
J_ bands are ^3 = ^ (0, cos6>, — sin6>, 1) and ^4 = 
(0, — cos 9, sin 9,1) , where we write k x = fc cos 9 and 
k y = k sin 9. The linear combinations "fi + "J 2 and $3 + "J 4 
have non-zero amplitude only on A and D sites respectively. 
Any other state will break up into fast (J + ) and slow (J_) 
fermionic excitations, analogous to fast and slow modes in a 
birefringent medium. 

Staggered potentials: Staggered on-site potentials are a nat- 
ural perturbation to Hk in the context of cold atoms on an op- 
tical lattice. We can write the most general form of such a 
potential as 

A = ^1 [ A ° /4 + Al ^° + A 2(*7 1 7 3 + «7 2 7 5 ) 
fc 



+A 3 (z7 1 7 3 -*7 2 7 5 )]V'fe, (12) 



where we may set Ao = since this just corresponds to a 
uniform shift of the chemical potential. The Ai term violates 
chiral symmetry in the usual way but is Lorentz invariant and 
hence introduces a gap in the dispersion of the fermions 



E k = ±jA 2 + AJ 2 ± k 2 . 



(13) 



When (3=1 there are fiat bands at E = ±Ai that intersect the 
J + bands only at (k x , k y ) = (0, 0). The birefringence prop- 
erty discussed above is unaffected by the Ai term. We com- 
bine Z7 1 7 3 and i7 2 7 5 into a Lorentz invariant term (A2) and 
a Lorentz violating term (A3). There are two cases in which 
we have obtained simple analytic solutions for the spectrum: 
case I): Ai ^ 0, A 2 ^ 0, A 3 = 0, for which 



E k = 



A 2 ± v /(A 1 +A 2 ) 2 +4j2fc2 
-A 2 ± J(A 1 -A 2 ) 2 +4J 2 k 2 



4 



and case II): Ai ^ 0, A 2 = 0, A 3 ± 0, for which 



E k 




A 3 ) 2 + 4J 2 fc 2 +4J 2 £; 2 



(A 1 + A 3 ) 2 



4J 2 fc 2 +4J 2 fc 2 



In case I) the dispersion remains isotropic in momentum space 
and there are fiat bands when (3 — 1, whereas in case II), the 
dispersion becomes anisotropic, with the anisotropy governed 
by /3 through J± . In both cases, there is a shift in the spectrum 
and there will be at least one set of massive modes (however 
in both cases there can be a set of massless modes whose dis- 
persion is given by the upper half of a cone if Ai = ±A 2i3 
and A = tA 2 , 3 ). 

Interactions: as we consider spinless fermions, there will 
be no on-site Hubbard interaction, so we consider nearest 
neighbour interactions of the extended Hubbard type (for cold 
atoms in an optical lattice these will generally be weak): 



Hint = ^2 VijUiUj 



(14) 



Setting all of the Vij = Vq, we can write the interaction 
Hamiltonian in terms of spinors as 



Vb 
16 



(15) 



with %pf* = i/'It - The identity and 7 that appear in the ker- 
nels of the quartic interaction terms are the only elements of 
the Clifford algebra that either commute or anticommute with 
all of the elements of the Lorentz group and the chiral group, 
ensuring that the interactions remain invariant under any ro- 
tation of the lattice by the Lorentz group or relabelling of the 
plaquette indices by the chiral group. At the mean field level, 
the ('tptp) 2 term breaks the chiral symmetry by introducing an 
effective mass term mo 7 °, and the (■0 7 °V') 2 term renormal- 
izes the chemical potential as SI4 and is otherwise uninterest- 
ing. In the limit of weak interactions, the mean field interac- 
tion Hamiltonian is: 



H^ = £ 4{ [(Sh + mo7°) + (m l7 V + m 2 / 7 2 
fe 

+m 3 i 7 3 + m 5 i 7 5 )] Vfc, (16) 

where S = (ha) + (tib) + (nc) + (no), and the order 
parameter for staggered charge density wave order mo = 
(tia) — { n B) — (nc) + { n D) arise from the Hartree term. The 
remaining masses, mi, m 2 , m 3 and arise from the Fock 
term - if these are dropped and (3 = 0, we recover the mean- 



field approximation of the Gross-Neveu model 111 611 . Similarly 
to a Ai 7 ° staggered potential, the Hartree term leads to mas- 
sive excitations, but does not destroy fermion birefringence. 
The detailed study of interactions when (3 ^ is a topic for 
future investigation. 



For small values of /3, when H m t is added to Eq. (01 there 
is a mapping between the weak interaction strength regime 
considered above to the strong interaction strength limit that 
preserves the property of birefringence: 



E k (p,V )=!3E k (p- 1 ,p- 1 Vo). 



(17) 



This arises from the appearance of the chiral symmetry gen- 
erators, 7 3 and 7 5 in the kinetic energy and their duality with 
Lorenz group generators 7 ° 7 1 , 7 ° 7 2 - Upon choosing a differ- 
ent representation of Clifford algebra elements, one can trans- 
form 7 ° 7 1 <-> 7 3 , 7 ° 7 2 o 7 5 . 

Topological defects: broken chiral symmetry at f3 ^ 
implies that there cannot be vortices, but domain walls 
of the form Ai(x) 7 ° where lim 2: _i. 00 Ai(x) = A and 
lim 2; _j._ 00 A±(x) = —A can occur. If (3 = then the form 
of the solutions with energy |e| < A is well known. When 
(3 we can find zero energy bound states with different 
spatial extents for the + and - solutions: 



where u + = (1, 0, i, 0), u_ = (0, —i, 0, 1), and n± = 1/2J±. 

In this Letter we have demonstrated a model whose low 
energy excitations are birefringent fermions that arise from 
broken chiral symmetry. We discuss the low energy proper- 
ties of the model and illustrate the meaning of broken chiral 
symmetry in our model. We argue that such a model could be 
realised by cold atoms in an optical lattice. This might not be 
the only route - as noted in a similar context in Ref. Hill an- 
other approach might be through an appropriately engineered 
semiconductor heterostructure. 

An important feature of the birefringent fermion dispersion 
that we find here is that the slopes of the J + and J_ bands can 
be controlled by the parameter (3. In particular when (3=1, 
there can be flat bands in the spectrum and these flat bands 
are robust to the addition of a staggered potential Ai 7 ° and 
weak nearest neighbour Hubbard interactions at the Hartree 
level. Flat bands such as Landau levels can lead to interest- 
ing correlated phases when interactions beyond mean field are 
taken into account lfl7ll . This suggests that future avenues for 
research on this model could include the study of such corre- 
lated phases when (3 ^ 0, and the generalization of the model 
to fermions with spin. Including spin would allow for on-site 
Hubbard interactions, which would be considerably compli- 
cate matters and require techniques similar to those that have 
been used to study high temperature superconductors 1 18]. 
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